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OF THE 
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ABSTRACT 

The albedo and dichroic polarization of spheroidal 
particles are studied. For the description of the phase func- 
tion two asymmetry parameters g^ | and g± characterizing 
the anisotropy in forward/backward and left/right direc- 
tions are introduced and calculated. The consideration is 
based on the solution to the light scattering problem by the 
Separation of Variables Method (Voshchinnikov and Fara- 
fonov, 1993). 



1. INTRODUCTION 

In many scientific and engineering applications pro- 
late and oblate spheroids are appropriate models for real 
particles. We consider the light scattering by homoge- 
neous spheroids using the Separation of Variables Method 
(S VM). The optical properties of prolate spheroids of var- 
ious aspect ratios a/b for several refractive indices m are 
calculated and the results for the particles of the same vol- 
ume are compared. 



2. GENERAL DEFINITIONS AND METHOD 

A spheroid (ellipsoid of revolution) is obtained by 
the rotation of an ellipse around its major axis (prolate 
spheroid) or its minor axis (oblate spheroid). The ratio 
of the major semiaxis a to the minor semiaxis b (i.e. the 
aspect ratio a/b) characterizes the particle shape which 
may vary from a nearly spherical one {a/b « 1) to a 
needle or a disk (a/b ^ 1). 

We assume that an incident plane wave has the 
wavelength A. Let a denote the angle between the 
propagation direction and the rotation axis of the spheroid 

(0° < a < 90°). 

For the axial propagation (a — 0°), there is no 
polarization of transmitted radiation due to symmetry. If 
a 7^ 0°, two cases of polarization of the incident radiation 
have to be considered: the electric vector E is parallel (TM 
mode) or perpendicular (TE mode) to the plane defined 
by the spheroid's rotation axis and the wave propagation 
vector. 



The size parameter is given by 
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where ry is the radius of the sphere whose volume is equal 
to that of the spheroid. The radius ry for prolate spheroids 
is defined as 
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One usually calculates the efficiency factors Q = 
C/G which are the ratio of the corresponding cross- 
sections C to the geometrical cross-section G of the prolate 
spheroid (the area of the particle's shadow) 



G{a) — nb (a^ sin^ a + b^ cos^ i 
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In order to compare the optical properties of the parti- 
cles of different shapes it is convenient to consider the ra- 
tios of the cross-sections for spheroids to the geometrical 
cross-sections of the equal volume spheres, C/irry. For a 
prolate spheroid they can be found as 
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The albedo of a particle can be calculated from the 
extinction and scattering cross-sections 
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In general, the radiation scattered by aligned 
spheroidal particles has an azimuthal asymmetry that pro- 
vokes a non-coincidence of the directions of the radia- 
tion pressure force and of the wave-vector of incident ra- 
diation (Voshchinnikov, 1990; Il'in and Voshchinnikov, 
1998). Another consequence of the azimuthal asymme- 
try is the anisotropy of the phase function in the left/right 
direction. The geometry of the phase function in for- 
ward/backward and left/right directions may be charac- 
terized by two asymmetry parameters g^ and g±, respec- 
tively. Expressions for them can be found from the consid- 
eration of radiation pressure (Voshchinnikov, 1990; Il'in 
and Voshchinnikov, 1998) 
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Figure 1: Normalized scattering cross-sections for elongated spheroidal particles in dependence on size parameter 27ra/A. 
Transformation to the parameter xv is given by the expression xy = {a/b)^'^'^{27ra/\) 



A{m,XY, a/b, a) sin a — B{m,XY, a/b, a) cos a 
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(7) 
where Qsca is the scattering efficiency factor, and the co- 
efficients A and B are 

A~ IC{xy,a/b,a) / / 1(6*, (^) cos0sin6'd0d93, 
Jo Jo 

(8) 

B = IC{xy,a/b,a) / / i(0, 1^9) sin'^ 0cos(^d0d(y9. 
Jo Jo 

(9) 
In Eqs. (H) - (||), /C is a parameter, i(0, 1^9) the dimension- 
less intensity of scattered radiation (phase function). From 
a symmetry consideration, it is clear that for spheroids 
gj^ = if a = 0° or 90°. Note that in all cases the follow- 
ing inequality is valid 



1 < 
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The dichroic polarization efficiency is defined by the 
extinction cross-sections for TM and TE modes 
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This ratio describes the efficiency to polarize light trans- 
mitted through an uniform slab consisting of non-rotating 
particles of the same orientation. 

The optical properties of spheroidal particles can be 
determined by various methods of light scattering theory 



(see Mishchenko et al., 2000 for a review). We use the 
SVM's solution developed by Farafonov and numerical 
code based on it (see Voshchinnikov and Farafonov, 1993 
for more details). A comparison of different numerical 
codes and benchmark results can be found in the paper of 
Voshchinnikov et al. (2000). 

The main problem of the SVM for spheroids is the 
difficulties with computations of the spheroidal wavefunc- 
tions. Especially it is related to very elongated particles 
with sizes larger than the wavelength because the standard 
expansions of the prolate spheroidal wavefunctions in se- 
ries of the Legendre functions do not converge. Farafonov 
and Voshchinnikov (2000, in preparation) have considered 
new expansion of the prolate wavefunctions that opens a 
possibility to calculate the optical properties of particles 
with a/b ^ 1. Firstly, this method was applied to calcula- 
tions of the radial wavefunctions with the index m = 1 that 
allows us to study the case of axial propagation of radia- 
tion. Some results are shown in Fig. |l]for non-absorbing 
spheroids with the refractive index m ~ 1.7 + G.Oi. As 
the extinction and scattering cross-sections were the same 
(with 6 and more digits), the values presented are expected 
to be correct. 

3. NUMERICAL RESULTS 

We present some results illustrating the behaviour of 
the optical properties of prolate spheroidal particles in a 
fixed orientation. 



3.1 Albedo 

The integral scattering properties of particles are char- 
acterized by their albedo. This quantity depends on the 
particle size and, in general, on the particle shape. Figure ^ 
shows the size dependence of the albedo for spheroids with 
m = 1.3 + 0.05i and m = 1.7 + 0.7i and the aspect ratios 
a/6 = 2 and 10. The calculations were made for pro- 
late particles and a — 0° (we adopt that the incident radi- 
ation is non-polarized). In this case in comparison with 
others (oblate spheroids, oblique incidence of radiation) 
the largest deviations of the ratio A(spheroid)/A(sphere) 
from unity occur (see Voshchinnikov et al., 2000 for dis- 
cussion). 





Figure 2: Albedo of spheroidal particles normaUzed rela- 
tive to albedo of spherical particles with the same refractive 
index 

It is seen that the albedo for large non-spherical par- 
ticles becomes close to that of spheres. Our calculations 
made for particles with different absorption show that 
the distinction of the albedo for spheres and spheroidal 
particles remains rather small (within ~ 20 %) if the ratio 
of the imaginary part of the refractive index to its real part 
k/ii k, 0.2 - 0.3. 

3.2 Asymmetry parameter 

Another characteristic of scattered radiation is the 
asymmetry parameter describing the spatial distribution 
of scattered radiation around a particle. Usually the 
anisotropy in forward/backward direction is only consid- 
ered. However, the radiation scattered by any aligned 
non-spherical particle possesses also an anisotropy in the 
right/left direction in the case of oblique incidence of radi- 
ation (see Fig. ^. 



As it is seen from Fig. ^both asymmetry factors g^ | and 
g± change with a/b and a. The values of radial asymmetry 
factor 5|| decrease with a growth of a when the path of 
radiation reduces from 2a (a — 0°) to 26 {a — 90°). The 
transversal asymmetry factor g± can be rather large and 
even exceeds the radial one. Because the geometry of light 



Figure 3: Geometry of light scattering by a prolate 
spheroid. The wave-vector of incident non-polarized ra- 
diation forms the angle a = 45° with the rotation axis of 
the spheroid. Short-dashed curve shows the angular distri- 
bution of the scattered radiation (phase function). The part 
of scattered radiation which is symmetric relative to the 
direction of incident radiation is plotted by long-dashed 
curve. The radial and transversal asymmetry factors and 
their values are indicated 



scattering by very elongated spheroids approaches that of 
infinite cylinders,^] such particles scatter more radiation 
"to the side" than in forward direction. 

The size dependence of asymmetry factors is plotted 
in Fig. g|. It shows that the variations of g||(xv) for 
spheroids with a/6 = 2 are rather similar to those of 
spheres. However, the radiation scattered by spheroids has 
a noticeable azimuthal dependence which is absent for 
spherical particles. If xy ^ 4 the azimuthal anisotropy of 
scattered radiation reduces and g± drops. 

3.3 Polarization 

If a volume contains aligned non-spherical particles, 
the initially non-polarized incident radiation will be 
partially polarized after having passed the volume. The 
simplest and at the same time extreme case of particles' 
alignment is the perfect alignment of non-rotating particles 
(picked fence orientation). The maximum polarization 
usually occurs when the major axes of the particles are 
perpendicular to the direction of the incident radiation 
{a = 90°). 

The behaviour of the polarization efficiency P/t for 
non-absorbing and absorbing spheroids is shown in Fig. H. 
It is clearly seen that a relatively large particles produce 
no polarization independent of their shape. For absorb- 
ing particles, it occurs at smaller xy values than for non- 
absorbing particles. This effect should depends on the 



' In this case the scattered radiation forms the conical surface with the 
opening angle 2a. 
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Figure 4: Angular dependence of the asymmetry factors 
of scattered radiation 



imaginary part of the refractive index as well. 



4. CONCLUSIONS 

The albedo of large non-spherical particles exhibits 
only a weak dependence on the particle shape if the ratio 
of the imaginary part of the refractive index to its real part 

fc/n ;> 0.2 - 0.3. 

The radiation scattered by aligned spheroidal particles 
has an azimuthal asymmetry and its geometry may be 
described by two asymmetry parameters g|| and g± 
showing the deviations from the symmetric scattering in 
forward/backward and left/right directions, respectively. 
The transversal asymmetry factor g± can be rather large 
and even exceeds the radial one, therefore, very elongated 
spheroids scatter more radiation "to the side" than in 
forward direction. 

Particles larger than a certain minimum size do not po- 
larize the transmitted radiation independent of their shape. 
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Figure 5: Size dependence of the asymmetry factors of 
scattered radiation 
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Figure 6: Size dependence of the polarization efficiency 
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